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The internal dynamics of active gels, both in artificial (in-vitro) model systems and inside the 
cytoskeleton of living cells, has been extensively studied by experiments of recent years. These 
dynamics are probed using tracer particles embedded in the network of biopolymers together with 
molecular motors, and distinct non-thermal behavior is observed. We present a theoretical model of 
the dynamics of a trapped active particle, which allows us to quantify the deviations from equilibrium 
behavior, using both analytic and numerical calculations. We map the different regimes of dynamics 
in this system, and highlight the different manifestations of activity: breakdown of the virial theorem 
and equipartition, different elasticity-dependent ’’effective temperatures” and distinct non-Gaussian 
distributions. Our results shed light on puzzling observations in active gel experiments, and provide 
physical interpretation of existing observations, as well as predictions for future studies. 

PACS numbers: 


I. INTRODUCTION 


In-vitro experiments have probed the non-thermal (ac¬ 
tive) fluctuations in an ’’active gel”, which is most com¬ 
monly realized as a network composed of cross-linked 
filaments (such as actin) and molecular motors (such 
as myosin-II) [IHl]- The fluctuations inside the ac¬ 
tive gel were measured using the tracking of individual 
tracer particles, and used to demonstrate the active (non¬ 
equilibrium) nature of these systems through the break¬ 
ing of the Fluctuation-Dissipation theorem (FDT) [2]. 
In these active gels, myosin-II molecular motors generate 
relative motion between the actin filaments, through con¬ 
sumption of ATP, and thus drive the athermal random 
motion of the probe particles dispersed throughout the 
network. This tracking technique was also implemented 
in living cells [3 H US]. The motion of these tracers 
in cells was also shown to deviate from simple thermal 
Brownian diffusion. 

There are several puzzling observations of the dynam¬ 
ics of the tracer particles inside the active gels, for ex¬ 
ample, the distinct non-Gaussianity of the displacement 
correlations and their time-dependence [BITIIHI- We pro¬ 
pose here a simple model for the random active motion of 
a tracer particle within a (linearly) elastic active gel, and 
use our model to resolve their distinct non-equilibrium 
dynamics. On long time-scales the tracer particles are 
observed to perform hopping-like diffusion, which is be¬ 
yond the regime of the present model, and will be treated 
in a following work, as will be the introduction of non¬ 
linear elasticity [3| . The activity is modeled through col¬ 
ored shot-noise unmi] and the elastic gel is described 
by a confining harmonic potential. We use the model to 
derive expressions directly related to the experimentally- 
accessible observations: such as the position and veloc¬ 
ity distributions and their deviations from the thermal 
Gaussian form. Our model allows us to offer a physi¬ 
cal interpretation to existing experiments, to character¬ 
ize the microscopic active processes in the active-gel, and 


to make specific predictions for future exploration of the 
limits of the active forces and elasticity. The simplicity of 
this model makes this model applicable to a wide range 
of systems, and allows us to gain analytic solutions, intu¬ 
ition and understanding of the dynamics, which is usu¬ 
ally lacking in out-of-equilibrium systems. This would be 
more difficult to obtain with more complex description of 
the gel, such as visco-elastic that has more intrinsic time- 
scales. 


II. MODEL 


Our model treats a particle in a harmonic potential, 
kicked randomly by thermal and active forces (active 
noise) m- The corresponding Langevin equation for the 
particle velocity v (in one dimension or one component 
in higher dimensions, with the mass set to m = 1) 


i) 


-Au + fa + fr 


dU{x) 

dx 


( 1 ) 


where A is the effective friction coefficient and the har¬ 
monic potential is: U{x) = kx"^/2, with k proportional 
to the bulk modulus of the gel (related to the gel density, 
cross-linker density and other structural factors). The 
thermal force /t is an uncorrelated Gaussian white noise: 
{frit)frit')) = 2XT6{t — t'), with T the ambient temper¬ 
ature, and Boltzmann’s constant set to /cb = 1. 

We model the active force fa as arising from the in¬ 
dependent action of molecular motors, each motor 
producing pulses of a given fixed force ±/o, for a dura¬ 
tion At (either a constant or drawn from a Poissonian 
process with an average value At, i.e. shot-noise), with 
a random direction (sign). The active pulses turn on ran¬ 
domly as a Poisson process with an average waiting time 
T (during which the active force is zero), which deter¬ 
mines the ’’duty-ratio” of the motor (the probability to 
be turned ”on”): pon = At/(t -|- At). 
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III. RESULTS: MEAN KINETIC AND 
POTENTIAL ENERGIES 


The mean-square velocity and position fluctuations of 
the trapped particle, essentially the mean kinetic (T„ = 
(n^)) and potential {T^ = k{x‘^)) energies, can be calcu¬ 
lated for the case of shot-no ise force correlations (details 
given in the Appendix, Eqs Ala|^ and Figs 3 5). Note 
that the mean (•) is over many realizations of the sys¬ 
tem, or over a long time. In the limit of vanishing trap¬ 
ping potential the position fluctuations (a:^) diverge, but 
the potential energy approaches a constant: T^lk^o —>■ 
/qAt/A. The kinetic energy approaches the constant 
value for a free particle m- ^ r,/(l-bAAr). 

We therefore find that the virial theorem is in general 
not satisfied in this active system, which in a harmonic 
potential gives Ty\f.q = Tx\eq, even in the limit of weak 
trapping. The virial theorem, and equipartition, breaks 
down due to the strong correlations between the particle 
position and the applied active force: In the limit of per¬ 
fect correlations, the particle is stationary at a; = ±a:o 
when the force is turned on (the stationary position in 
the trap where the potential balances the active force: 
xq = fo/k), and at x = 0 when it is off. In this extreme 
case the potential energy is hnite while the kinetic energy 
is zero. 

In the limit of strong trapping k —>■ oo,fc/ A^ 1 , 
the potential energy behaves as: oc k~^ (Eqs |B6|B9| ), 

while the kinetic energy decays faster as: T„ cx k~'^''^ 
(Eq |D2| ). One can understand this limit as follows: When 
the trapping is very strong, the shortest time-scale in the 
problem is the natural oscillation frequency in the trap, 
Wfc ~ '/k. In this regime of fcAr^ ^ 1 we find that during 
the active pulse At, the particle reaches Xg, and the mean 
potential energy is therefore proportional to ~ kx^ oc 
1/fc. The kinetic energy in this limit decays faster, since 
the fraction of time that the particle is moving is only 
during the acceleration phase determined by the time- 
scale ~ y/k. We therefore find that in the presence of 
strong elastic restoring forces the potential energy will be 
much larger than the kinetic energy, in an active system 
(Tj, ^ Ty). This was recently found in the study of active 
semi-flexible polymers [T^ . 

Note that in a real active gel the different parameters 
maybe coupled: larger local density of the network fil¬ 
aments increases the local value of the elastic stiffness 
parameter fc, but may also increase locally the density 
of motors and their ability to exert an effective force, 
thereby increasing and /g. The tracer bead behavior 
as expressed by T„ and can therefore be a complex 
function of the local network parameters. 


IV. RESULTS: VELOCITY AND POSITION 
DISTRIBUTIONS 

The distributions of the velocity and position in the 
different regimes are shown in Figj^ for the case of a sin¬ 


Log(P(v)) X(x) Log(P(x)) 



FIG. 1: Distribution of position and velocity for particle 
trapped within a harmonic trap of various stiffness {k = 
1,1000 for a,b and c,d respectively), and for different regimes 
of friction (A = 50,0.1 for a,c and b,d respectively). The 
time-scale of the active bursts is At = 0.1, the amplitude of 
the active force /o = 1 and the waiting time t = 1 (so that 
Pon ~ 0.1). For simplicity we plot the behavior for the case 
of a single motor with a constant burst duration. The insets 
compare the simulated distribution (blue line) to the analytic 
approximation (red dashed line), in log-linear scale, as simple 
Gaussians or as a sum of shifted thermal Gaussians. 


gle active motor. The simulations of the model were car¬ 
ried out using explicit Euler integration of Eqj^ (see also 
the Appendix for details). We study this case in order 
to highlight the deviations from Gaussian (equilibrium¬ 
like) behavior, which is restored by many simultaneous 
motors m- In an infinite gel, with a constant density of 
motors, we may therefore treat the distant (and numer¬ 
ous) motors as giving rise to an additional thermal-like 
contribution to the tracer dynamics (Eqs D1|D2 ), while 
the non-equilibrium behavior is dominated by a single 
proximal motor [1]. 

In the limit of weak damping, A At <C 1, both the posi¬ 
tion and velocity distributions are very close to Gaussian, 
with the width of the Gaussia n distributions given by T„ 
and Tx (Figj^,d, Eqs B4|cT ): P{v) oc exp (—u^/(2T„)), 
P(x) oc exp b^fcx^/(22^Jb 


In the highly damped limit, A At 3> 1, the distribu¬ 
tions become highly non-Gaussian (Figj^,c). We can 
make a useful approximation in this limit, by neglecting 
the inertial term in Eq. Q and get the following equation 
for the particle position x inside the potential well 


\v = -kx + fa + It (2) 

^ \ , fa + It 

=> X = —Ax H-r- (3) 

A 

where A = fc/A. This equation is now analogous to the 
equation for the velocity x of a free particle (Eq. Q when 
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U{x) = 0). Due to this analogy we can use the analytic 
solutions for the free particle m to describe the parti¬ 
cle position in the well. For weak trapping (Figj^), we 
therefore expect the position distribution to be roughjj 
Gaussian, since we are in the limit of AAr ^ 1 of Eq.([^, 
with a width given by (from Eqj^S6) 


T' 


T, 


PonNmX (^XAt + e 
k'^Ar 

r, PonNm{AT) ^2 

^ 7 : V •'0 

A(l + A(Ar)j 


/o 


( 4 ) 

( 5 ) 


where describes the case of a constant At and the 
case of a Poissonian burst distribution, and fits well the 
calculated distribution (inset of Fig[^). In the limit of 
weak confinement we expect the velocity distribution to 
approach the behavior of the free damped particle 
which is well approximated as a sum of thermal Gaus- 
sians, centered at v = 0,±no (vq = Fg/X). This is indeed 
a good approximation, as shown in the inset of Figj^. 

For strong potentials (AAr 1, Figj^) we expect 
from the analogy given in Eq. @ that the spatial distribu¬ 
tion is now well described by the sum of shifted thermal 
Gaussians (Fi gS^)[ii] , centered at a: = 0, ixo- The ve¬ 
locity distribution in this regime is also non-Gaussian: 
the maximal active velocity is of order uq at the ori¬ 
gin of the potential, but since the particle immediately 
slows due to the confinement (up to a complete stop at 
±a:o), the peaks of the distribution are located at roughly 

±Vo/2. 


V. RESULTS: NON-GAUSSIANITY OF THE 
DISPLACEMENT DISTRIBUTION 


The distribution of relative particle displacements 
(Van Hove correlation function) P{Ax{ti^)), where 
Ax{Ta,) = x{t + To,) — x(t) (juj is the lag-time duration), 
is a useful measure for the particle dynamics. We plot 
it in Figj^ for the interesting regime of strong confine¬ 
ment and damping, and compared to the distribution of 
particle positions P[x) (Fi sS^)- We see that P{Ax{t^)) 
has double the number of peaks of P{x), and is distinctly 
non-Gaussian for all r^. In Fig|^ we show that the same 
qualitative behavior is obtained for Poissonian burst du¬ 
ration. 

The deviations from Gaussianity are quantified in 
Figj^ using the Non-Gaussianity Parameter (NGP) of 
the displacement distributions: k = {Ax"^)/3{Ax‘^)'^ — 1. 
This deviation of the kurtosis from the value for a Gaus¬ 
sian is an established measure for studying distributions 
P5] . We find that the NGP has a finite value for —>• 0. 
This is as a consequence of the periods during which the 
particle is accelerated by the active force, and the re¬ 
sult is a finite probability for displacements of the order 
of Ax ~ vqTi^ (inset of Figj^). With increasing Tuj the 



FIG. 2: (a) Distribution of particle displacements P{Ax{tu,)), 
for various lag time duration (blue lines), for a single mo¬ 
tor and constant burst duration. The traces correspond to 
increasing time lag durations (black arrow), in the range: 
1 > Tij > 5 X 10““^. The red dashed line denotes the spa¬ 
tial distribution P{x), and the black dashed lines denotes 
P{Ax{(x)) (Eq|^ for the approximation of P{x) as a sum 
of three Gaussians. Inset shows the displacement distribu¬ 
tion for very short lag times Tu,. Parameters as in FigdJ:- 

(b) Calculated NGP for the P{Ax{tui)), for various num¬ 

ber of motors {Nm), and conhnement strength. The short 
horizontal black lines denote: (left) the NGP of P{x), and 
(right) of P{Ax{oo)) (Eq|^, for the k — 1000, = 1 case. 

(c) Displacement distributions (as in (a)) for a calculation 
without the inertia term (Eqj^ using: k — 1000, Poissonian 
(At) = 0.1, and increasing lag time indicated by the arrow 
Tuj = 10“^, 2.5 X 10“®,5 X 10“®, 10“2,1), and (d) the corre¬ 
sponding NGP, comparing the simulation (solid grey lines) to 
the analytical result (see Appendix for details, dashed lines), 
for k = 300,1000 (top, bottom). The NGP with inertia is 
given by the solid black lines. 


NGP reaches a maximum, at lag times that are of order 
At, where the full effect of the active bursts is observed. 

We find that the maximal value of the NGP for 
P{Ax{t^)) is close to the NGP of P{x) (Figj^), which 
is a function for which we have a good analytic approx¬ 
imation [H] (Eqs |Fl|F2D . In the limit of —>■ oo the 
calculated NGP remains finite, and can be calculated an¬ 
alytically since the displacement distribution becomes 


/ OO 

P{x)P{x -\- Ax)dx 

-OO 


( 6 ) 


and the P(x) in this regime is well approximated by the 
sum of shifted thermal Gaussians (inset of Fig ^). This 
calculation fits well the simulated result (Fig^). For 
a larger number of motors, the distribution PJAx{ti^)) 
approach a Gaussian (Fig[^,S4). 

In the limit where we discard inertia from the equations 
of motion (Eq|^, we can calculate the NGP analytically 
(see details in Appendix). In Figj^ we plot the displace¬ 
ment distributions for this case, and in Fig[^ we show 
that indeed the analytical calculation describes exactly 
the simulation results. We find that this treatment cap- 
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tures correctly the qualitative features of the full system, 
such as the position of the peak, followed by a constant 
value at long lag times. The large discrepancy is in the 
limit of —>■ 0, where the inertial effects of the oscilla¬ 
tions inside the trap are missing from Eqj^ 


VI. RESULTS: FDT 


An alternative method to characterize the non¬ 
equilibrium dynamics is through the deviations from the 
FDT [ 2 ] . We can quantify these deviations by defining an 
effective temperature, using the Fourier-transform of the 
position fluctuations (S^xi^)) and linear response (sus¬ 
ceptibility of the position to an external force x(‘^)) of the 
system. We can calculate both for our trapped particle 
position using Fq.(|^ for the A At 1 limit, to get (for 
Poissonian burst duration At, see details in Appendix, 


Fqs |Gll[G3 ) 


Tpdt{^) 


_ NjjiPonfo (At) 

2Im(x(a.)) “ A(1 + (^(At))^) + ^ ^ 


Note that is independent of the shape of the 

harmonic potential (fc), and is identical to the result for 
a free active particle m- This result highlights the fact 
that while different ” effective temperatures” in an active 
system {Ty and Tx, Eqs|4|^ give a measure of the activ¬ 
ity, they can have very different properties. 


VII. DISCUSSION 


We now use our results to interpret several experiments 
on active gels in-vitro, and extract the values that charac¬ 
terize these active systems. In [2] the break-down o f the 
FDT was measured. Comparing to our Tp dt (Eq G3) 


we find that the onset of the deviation from equilib¬ 
rium occurs for frequencies w < At~^, from which we 
find that: At « 100 msec, which is the scale of the re¬ 
lease time of the myosin-II-induced stress [2] in this sys¬ 
tem. The measured deviation from the FDT was found 
to increase with decreasing frequency [2], and at the 
lowest measured frequencies the ratio was found to be 
Tpdt{'jJ —> 0)/T « 20 — 100. This number fixes for us 


the combination of the parameters given in Eq.(G3). 


Recent experiments shed more detail on the active mo¬ 
tion in this system [T], and it was found that the tracer 
particle performs random confined motion interspersed 
by periods of large excursions. The confined motion 
part can be directly related to the mean-square displace¬ 
ment in our model Tx (Eq[^, and is observed to be a 
factor of Tx/T « 10 — 50 larger than in the inert sys¬ 
tem (not containing myosins) These values are in 
general agreement with the values extracted above for 
TpDT from [ 2 ], and note that we predict (Eqs 5|G3 |: 
TporipJ —t 0)/Tx = 1-1- A(At) > 1. 

Furthermore, in these experiments [T] it was observed 
that the distribution of relative particle displacements 


P{Ax{tix)) is highly non-Gaussian. Comparing to Figj^ 
we note that both the experiments and in our calculations 
the NGP has a finite value for 0. With increasing 
the NGP reaches a maximum, both in the experi¬ 
ments and in our calculations (Figj^,d). By comparing 
to our model we expect the peak to appear at t^, « At, 
so the observations [T] suggest the burst duration is of 
order At « 1 — lOsec, in agreement with similar studies 
[ail]. Note that very similar NGP time-scales were ob¬ 
served in living cells [mill] Our model predicts that the 
maximal value of the NGP is a non-monotonous function 
of Pon, and this may be explored by varying the concen¬ 
tration of ATP in the system. Furthermore, from our 
model we predict that the NGP decrease with decreas¬ 
ing active force, and increasing stiffness of the confining 
network (Figj^,d, Eq F2). These predictions can be re¬ 
lated to the observed activity-dependence of the NGP in 
cells [H], and the decay of the NGP during the aging and 
coarsening of an active gel |7]. 

The large observed deviations from Gaussianity indi¬ 
cate that the particle is in the strong confinement regime: 
A (At) > 1. The maximal value of the observed NGP 
« 2 — 4 can be used to get an esti mat e of Tx, by tak¬ 
ing it to be equal to NGP^ax (EqF2). This gives us: 
Tx ^ 10 — SOksT and Pon « 2 — 3%. This value of Tx is 
in good agreement with the estimate made above. The 
value of Pon is in agreement with the observation that 
the waiting-time between bursts is much longer than the 
burst duration |2], and with the measured duty-ratio of 
myosin-II Ca¬ 
in the limit of t^^ —>■ oo the observed NGP of the 
displacement distribution P{Ax{tuj)) decays to zero [T], 
while for the calculated confined particle the NGP re¬ 
mains finite (Figj^,d). At long times (^ lOsec) the 
observed trajectory has large excursions [Tj, which we 
interpret as the escape of the particle from the confining 
potential. The ensuing hopping-type diffusion, causes 
the NGP to vanish, as for free diffusion [T^]. Within our 
model we therefore interpret the observed time-scale of 
the vanishing of the NGP, t^^ « 10— 100 sec, as the time- 
scale which corresponds to the mean trapping time of the 
bead within the confining actin gel. Beyond this time- 
scale the bead has a large chance to escape the conhne- 
ment, and hop to a new trapping site, which corresponds 
to a re-organization of the actin network. The real actin- 
myosin gel undergoes irreversible processes that make 
its properties time-dependent and render it inhomoge¬ 
neous ladMS]. Such effects make the comparison to 
the model much more challenging. Large deviations from 
Gaussianity were also observed for the Van Hove corre¬ 
lations in other forms of active gels [20]. 


VIII. CONCLUSION 

We investigated here the dynamics of a trapped active 
particle, with several interesting results: (i) The activ¬ 
ity leads to strong deviations from equilibrium, such as 
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the break-down of the virial theorem and equipartition. 
We find that in the presence of elastic restoring forces 
the activity is mostly ’’stored” in the potential energy of 
the system, (ii) Different ’’effective temperatures” give a 
measure of the activity, and some are dependent on the 
stiffness of the elastic confinement, (iii) The displace¬ 
ment, position and velocity distributions of the parti¬ 
cle are highly non-Gaussian in the regime of strong elas¬ 
tic confinement and small number of dominant motors. 
These distributions can be used, together with our sim¬ 
ple model, to extract information about the microscopic 
properties of the active motors. Note that in our model 
the activity affects the motion and position distributions 
of the trapped particle, which is complimentary to mod¬ 
els where the activity drives only the large-scale reorga¬ 
nization that moves the particle between trapping sites 


El 127], or leads to network collapse El- The results of 
this model are in good agreement with observations of the 
dynamics of tracer beads inside active gels, and the sim¬ 
plicity of the model may make it applicable for a wide 
range of systems. More complex visco-elastic relations 
can be used in place of the simple elasticity presented 
here, to describe the dynamics inside living cells E3^, 
as well as non-linear elasticity El- Note that in most cur¬ 
rent experiments on actin-myosin gels, the myosin-driven 
activity is strong enough to lead to large-scale reorgani¬ 
zation of the actin network, eventually leading to the 
network collapse EmZHISI- In order to observe the ac¬ 
tive motion for the elastically-trapped tracer in the intact 
network, which we have calculated, much weaker active 
forces will be needed. Our work can therefore give moti¬ 
vation for such future studied. 


Appendix A: Numerical simulations 


The simulations of the dynamics of the particle inside the ID harmonic potential were carried out using explicit 
Euler integration of Eq[2 We were careful to use a small time-step At, such that it was always an ord er of magnitude 
smaller than the smallest time-scale in the problem. The time-scales in the problem are: r,Ar and \/2jk, where the 
last time-scale is that of the oscillation frequency of the particle inside the harmonic potential. 

The iterative equations take the following form in terms of the sampling time At 

v{t + At) = v{t) + [—Au(t) — kx{t) + fa{t)] At + V2XTAtri, (Ala) 

x{t + At) = x{t) + v{t)At, (Alb) 


where 77 is a random Gaussian variable with zero mean and variance 1. Gonsidering that both the waiting time and 
the persistence time are exponentially distributed with mean values r and At, respectively, the iterative equation for 
the active force fa obeys 


fa{t + At) 


fait) 

if 

fait)^0 

prob. l — At/Ar 

fait) 

if 

fait)=0 

prob. 1 — At/r , 

0 

if 

fait)^0 

prob. At/A t , 

H-fJ} 

if 

fait)=0 

prob. At/r , 


where = {/,—/} with same probability. 


(A2) 


Appendix B: Position fluctuations of a trapped particle 


From the model equations of motion (EqQ, we can calculate the mean-square fluctuations in the particle position 
for a shot-noise force correlations with average burst duration At. We begin by Fourier transforming EqJ^to get 


— oj'^x = iojXx + fa + It — kx 


where the denotes the FT. From Eq B1 we get 


^ /a(w)+/T(w) 

^ = - 2 - ■ X I 7 

—— lUJA + k 


The fluctuations (correlations) are therefore 

{x^){uj) = {x{u!)x*{uj)) = 




(Bl) 


(B2) 


(B3) 
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FIG. 3: Calculated mean-square position fluctuatio ns (p lotted as a mean potential energy) for the trapped particle; Brown line - 
full solution, purple line- approximate solution (Eq B8l, blue line- approximate expression for the limit AAr ^ 1 (Eq|B9ll. 
In both panels we used Ar = 1, and (a) A = 10, (b) A = 0.01. In (a) the blue line agrees perfectly with the full solution, while 
in (b) it has a discrepancy at intermediate confinements. 


where we have: (/a)(w) = A^mPon/o i+(^at)^ (Poissonian shot-noise with mean burst length Ar), and (/|.)(w) = 2AT 
(thermal white noise) [llj . 

For the active part alone, we get 


^) = 


A^mPon/o 


1 


At 


27r Jg — kY + (wA)^ 1 -I- (wAt)' 


-dio 


(B4) 


The solution for this integral is quite lengthy. In the limit of weak trapping, fc —> 0, we get that the mean-square 
displacement diverges 


{x^) 


NmPonfo^T 


k {kAr -b A) 

such that the mean potential energy in this limit approaches a constant value 

T' ~ k{xY ^ 

' (fcAr -b A) 


(B5) 


(B6) 


In the limit of large k, we can expand the integrand of Eq.(B4) in powers of fc ^ to get this integral 

At 


ON n r'^ 

^ 2 \ _ mPonJ 0 / 


{xY = 


Iq fc2(l -b (wAt)^) 


du> 


(B7) 


which is also bound with a maximal frequency corresponding to the natural frequency of the harmonic trap. This 
integral gives a simple expression, which gives a good fit description as long as k ^ (Fig§ 


Y)k = arctan [^/k/2A^ 

TT rC 


(B8) 


Finding the value of k for whic h the scaling changes from {x^ k to ~ fc is simply by equating the large 
and small k limits of {x'^)k (Eq B8). 

In the limit of AAt ;g> 1 we find the simple approximate expression (Figj^ 


T' cs k{x^) cs 


ArN^nPonfa 
8A2(/cAT/2A-b 1) 


(B9) 


The numerical simulations, in the highly damped limit (AAt ^ I) indicate the k ^ and k 2 limits (Figj^). 
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FIG. 4: Calculated mean-square velocity fluctuations for the trapped particle: Brown line- full solution, purple line- approximate 
solution (Eq C31, blue line- highly damped limit (Eq |C4[ ), and the dashed blue line is the free-particle value m- In both panels 
we used At = 1, and (a) A = 10, (b) A = 0.01. 


Appendix C: Velocity fluctuations of a trapped particle 


Similar to the procedure for the position fluctuations described above, we ca n ca lculate the velocity fluctuations. 
The mean-square fluctuations in the particle velocity are given simply from Eq. (B4) by 


NmPonfo At 

2t 7o - fc)2 + (wA)2 1 + (wAt)2 


(Cl) 


The solution for this integral is again quite lengthy. As for the position distribution, we can find an approximation 
for the large k limit, using 


{v^) 


‘^^mPonfo 

TT 


r'^ uj^At 
'o k^l+iLoAT)^^^^ 


(C2) 


which is also bound with a maximal frequency corresponding to the natural frequency of the harmonic trap. This 
integral gives a simple expression, which gives a good fit description as long as fc 3> 

~ arctan [v^^Ar]^ (C3) 

The scaling of {v'^)k changes from (u^) ^ to {v"^) ~ k~^l'^ as k increases (Fig|^). 

In the limit of AAt ^ 1 we have the simple approximate expression 


ArN^-Ponfl 

4(A(1 -f AtA) -I- ttAt^ yjk^!%) 


(C4) 


which fits quite well the full expression in Fig|^. 

The numerical simulations, in the highly damped limit (AAr ^ 1) indicate the and limits (Figj^). 


Appendix D: Effective temperature due to forces from distant (and numerous) motors 


In a linear elastic medium, the displacements and stresses decay from a point source (at least) as 1/r^. Since 
there are numerous distant motors affecting the bead, their cumulative random forces are most likely to give rise 
to Gaussian distribution of position and velocities for the trapped particle. Each shell (of thickness dr) at radius r 
from the tracer beads has Nm{r) = Airr'^pdr motors (at constant density p), and therefore they contribute to the 
mean-square velocity the following contribution (in the limit of AAr 3> 1, using Eq C4) 


~ Nmir) ( NmPonfo^ 


At 


1 


4(A(1-b AtA)- b ttAt^v^) r^ 


(Dl) 






















FIG. 5: Simulated mean-square particle displacements (a) and velocity (b) in the limit of AAr ^ 1, using: At = 1, A = 50, 
fo = 1, Pan = 0.1. The dashed lines indicate the power-laws with exponents —1, —2 and —3/2 in (a) and (b) respectively. 





Ax(t„) 



FIG. 6: Simulated particle position distribution P{x) (red lines) and displacement distributions P{Ax{tu,)) (blue lines), for 
increasing number of motors: Nm = 5,10, 20 (left to right), using: Ar = 1, A = 50, fo = 1, Pon = 0.1. 


where we isolated the number of motors and the r-dependence of the active forces, and introduced a length-scale a 
beyond which the far-held calculation holds. Integrating this expression we get 

{v‘^)far ^ {v'^)o (dTrpa^) (D2) 

where (v^}o is the value for the single proximal motor given in Eq |C4[ We hnd that the far-held contribution of the 
distant motors is proportional to their density p. 


Appendix E: Displacement distribution for numerous motors 

As the number of motors kicking the particle {N„i) increases, we hnd that the distribution of the particle position 
becomes more Gaussian, even in the limit of larger damping A At ^ 1 and strong conhnement AAr 3> 1. We demon¬ 
strate this in Figj^ which shows that the position distributions P{x) and the displacement distributions P{Ax{tui)) 
approach a Gauassian for Nm larger than ~ 10. 


Appendix F: NGP for the highly damped limit 


We hnd that the maximal value of the NGP for P{Ax{t^)) is close to the NGP of P{x) (Figj^), which is a function 
for which we have a good analytic approximation m, given by (for a single motor) 


4(1 - 3po„)pl„Tl^ 


NGP{Nm = 1) 


3(1 + 2p2„r,.i)2 


(FI) 
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where is the effective temperature of the spatial distribution (Eq|^ for pon = 1- The maximal value of the 
NGP for a single motor, as a function of Pon is obtained from Eq|FT| at: Pon = a/(2 + 6a), and is given by 


NGPn 


(1 + 3a)^ 
3a(2 -t- 3a) 


(F2) 


where a = ksT/T^p. This is a monotonously decreasing function of the stiffness fc, due to the decrease in Txp in 
stiffer gels (Eqs.4,5). 


Appendix G: Effective temperature from the FDT, Tfdt 

Following m, and using EqU we can write for the AAr ^ 1 limit (when T = 0) 

Response : Xxx(oj) = 

Fluctuations : Sxx{^) = 

^ Tfdt(w) 


1 

7(fa; — A) 
fo 


(At 


A(A2 1 + (a;(AT))2 

_ NmPonfoi^T) 

2Im(x(a;)) A(1 + (a;(AT))2) 


resulting in Eq G3 


(Gl) 

(G2) 

(G3) 


Appendix H: Analytic calculation of the NGP without inertia 


To compute the expression of the NGP, we derive the mean quartic displacement (MQD) 
where it is time translational invariant 


^Ax^) in the regime 


{Ax^) = {A4) + {Axi)+6{Axl){Axl) , 


(HI) 


where the subscripts T and A refer respectively to the thermal and active contributions. The expression of the MSD 
is given by 


(A4) (t) 


2k^T 

k 



(Axl) (t) 


2kBT/^/k 
{r/TrY - 1 


(l - e-*/^) + e-*/^' - 1 


(H2) 


where Ti- = A/fc is a thermal relaxation time scale. The MQD under purely thermal conditions is related to the thermal 
MSD since the thermal process is Gaussian (Ax^) = 3 (Ax^) . To compute the active MQD, we separate the position 
displacement Ax^Yijtf) = x^^tf) — Xj^Yi) in several contributions, such that (Ax^) is a power law combination of 
these contributions. We compute each term using the active force statistics, and take the limit of large ti at fixed t 
corresponding to the time translational regime. The advantage of the separation we propose is that each term of the 
active MQD converges in such limit. The appropriate separation is 


AxA.,a (ti, tf) 


AxA,b (ti 5 tf) 


-^) J dt'x(ti - t')/A(t') , 

t 

J dt'x(t-t')/A(ti+tO , 


(H3a) 

(H3b) 
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where x(t) = e is the non-causal response function, and t = t{ — U is the time lag. In the time translational 

regime, we compute 


(t) = Tl 
(Aa^i^aAxA.b) (t) = Tl 
(t) = Tl 


3rr‘(2To + rr)(To + r)e (e*Ar _ i) 

A2(rr + T)(rr + 3T)(Tr(To + t) + 2tqt) 

3r^^r(2ro + Tr)(ro +T)e“^ (e*Ar — 

A2(Tr - T)(Tr + t){t^ + 3t){t,:{to + t) + 2tot) 

T/ (e Ar _ 1) e -r ^ 


(H4a) 

(H4b) 


A2(rr - ro)(r - Tr)(rr + r)2(rr(ro + r) - 2ror)(rr(ro + t) + 2tot) 

X ^4ro (t - Tr){T, + T)e^~^ + (r^ - to)(t - Tr) (tKto + rf - 4 tIt^) + 
+(t-o - Tr)(To + r)(Tr + r)e‘A (4r2r - r2(ro + r)) 

-2(to + r)e*A<- (2Torr + r(ro - rr)(2To + Tr)) (2tot - Tr(ro + r))^ , 


(Aa^A.aAxi^b) (0 = -Tj 


Srle-^ - l) 


T^(2to - Tr)(To + t) 


A2T(rr + r) V (Tr - - T)(ro(rr - 2t) + TrT) 

2TQe Tr + x) 


(t-Q - 'rr)(T-0 + Tr)(Tr(ro + t) - 2 tot ) 
(to + t)(t(to + Tr) - rorr)e^“- 


(Tr - To) (t2 - t2) 


(Axi^o) (0 = Tl 


(to + Tr) (r2 - 4rrr + 3 t2) 

3T,^e“^ / r(2ro - rr)(To + T)(Tr + t) 


A2T(Tr + t) V (Tr - 3r)(rr - T)(ro(rr - 2t) + TrT) 

8to (Tr + 

(Tr - To)(ro + Tr)(Tr(ro + t) - 2roT)(rr(ro + t) + 2Tor) 

4(to + r)(T(To +Tr) - roTr)e^“T _ .^^.^)gt(---) 


(tq + Tr) (r^ - 4rrT + 3 t2) 

T(2to + Tr){To + Tje'^r ^ 2Te^>- 


(to - Tr)(rr - r)(rr + 3r) 


(Tr + 3T)(Tr(ro + t) + 2roT) T -TrJ 


(H4c) 


(H4d) 


(H4e) 


from which we deduce 


(Axi) = (Axi_a) + 3 (Aa;l_aAxA.b) + 6 (Axl^^Axl^^) + 3 (A^A.aAxi^o) + (Axi_b) 


(H5) 
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